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Abstract 

In this paper we describe a strategy to study the Anderson model of 
an electron in a random potential at weak coupling by a renormaliza- 
tion group analysis. There is an interesting technical analogy between 
this problem and the theory of random matrices. In d — 2 the random 
matrices which appear are approximately of the free type well known to 
physicists and mathematicians, and their asymptotic eigenvalue distri- 
bution is therefore simply Wigner's law. However in d = 3 the natural 
random matrices that appear have non-trivial constraints of a geometri- 
cal origin. It would be interesting to develop a general theory of these 
constrained random matrices, which presumably play an interesting role 
for many non-integrable problems related to diffusion. We present a first 
step in this direction, namely a rigorous bound on the tail of the eigen- 
value distribution of such objects based on large deviation and graphical 
estimates. This bound allows to prove regularity and decay properties 
of the averaged Green's functions and the density of states for a three 
dimensional model with a thin conducting band and an energy close to 
the border of the band, for sufficiently small coupling constant. 



This paper is dedicated to Claude Itzykson. 



1 Introduction and Warning! 

The reader might be slightly surprised (perhaps even baffled?) by this con- 
tribution. Indeed it has nothing to do directly with the main topic of these 
Proceedings, namely strings, membranes and duality. It is however a tradition 
(and a good one!) to include in this cycle of conferences some applications of 
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quantum field theory to condensed matter physics. But in the case of our con- 
tribution, not only its subject, the mathematical analysis of disordered systems, 
is quite marginal from the point of view of the main theme of the conference, 
but also its form. It is not a "pedagogical introductory review" on the subject 
(for such a review we refer e.g. to [A][S]). We chosed indeed to give here the 
main technical elements of the proof of a non-perturbative regularity theorem 
on the Anderson model of an electron scattered by a random potential. We 
hope that glancing at this paper might nevertheless be amusing or even use- 
ful to the reader for two reasons. First the arid details of the mathematical 
bounds below will remind us how much effort is necessary to transform into 
proof even what certainly looks as the most trivial of all theoretical issues. Sec- 
ond we hope that the reader might be interested by the tool that was recently 
introduced in this area, namely a new category of random matrices which have 
constraints of a geometrical origin. The free random matrices with independent 
entries considered by Wigner, Dyson and others (and to the theory of which 
Claude Itzykson contributed so successfully) held the key to many unexpected 
theoretical developments which range from nuclear physics and confinement to 
quantum gravity and non-perturbative strings. There is now a general mathe- 
matical frame for such matrices, namely the free non commutative probability 
theory of Voiculescu [V]. In this frame, Wigner's law appears very general: it 
is simply the non commutative analog of the Gaussian law of large numbers for 
independent ordinary (commutative) random processes. But the theory of con- 
strained non commutative processes remains to be developed. There are many 
three dimensional problems (the BCS theory, the Anderson model, scattering 
of particles in 3 or 4 dimensional Minkovski space) that can be discretized as 
problems of random matrices with geometrical constraints of the type consid- 
ered below in section 4. These matrices seem definitely to have quite different 
scaling properties from the usual ones with independent entries. Perhaps they 
might also become useful tools in many other domains, when integrability and 
two dimensional conservation laws are no longer available. Therefore we think 
that the scaling properties of such random matrices of non-free type deserve to 
be further investigated, even numerically. 

2 The Anderson Model 

The Anderson model of a single electron in a random potential was introduced 
to study the concept of localization. 

Rigorous results in dimensions not restricted to 1 up to now concern the 
localization regime. The first result of this kind was the multiscale analysis 
of Frohlich and Spencer, who proved absence of diffusion for large disorder or 
energy out of the conduction band [FS] . It was later proved that in this regime 
localization holds, namely the expectation value of the modulus of the Green's 
function decays exponentially and the spectrum of the Hamiltonian is almost 



2 



surely pure point [DLS] , [FMSS] . The proofs have been simplified more recently 
by an argument due to Aizcnman and Molchanov [AM], which is based on 
considering the expectation value of the modulus of the Green's function to a 
power strictly smaller than one. 

However almost nothing rigorous is known about the "extended states" 
regime where the energy lies in the middle of the spectrum of the free Hamil- 
tonian and the disorder is weak. This is the regime where diffusive conduction 
should take place, at least in 3 dimensions. There are only results at d = oo by 
Klein [K]. 

We think that the program of renormalization group around singular man- 
ifolds should be useful to study this regime. This program was originally de- 
veloped for the study of many Fcrmions and in particular for a rigorous un- 
derstanding of the BCS theory of superconductivity [FT], [FMRT1]. In that 
context, the singular manifold was the Fermi surface. In the case of the An- 
derson model, the singular manifold is simply defined by the equation p 2 = E, 
namely it is a sphere for a free rotation invariant Hamiltonian; in the case of 
a lattice model breaking rotation invariance, one would have a slightly more 
complicated manifold such as e.g. 5Z i=1 (2 — 2cospi) = E. 

This renormalization group analysis consists in slicing the free Green's func- 
tion around the singular locus in momentum space, and performing a phase 
space analysis (by means of cluster expansions) which allows to compute the 
long range behavior governed by the effective theory near the singularity. 

In order to control this renormalization group analysis it is useful to un- 
derstand the behavior of the theory in a single momentum slice, hence slightly 
away from the singularity, and to prove that perturbation theory in this case is 
meaningful (i.e. asymptotic), using the fact that the coupling is weak to control 
the cluster expansion. 

In a previous paper [P] this single slice analysis was performed in two dimen- 
sions. In that situation the main problem consists in perturbing the identity 
matrix by a random matrix with independent coefficients, and this problem can 
be controlled using e.g. the asymptotics to the Wigner's law for the density of 
eigenvalues of this kind of matrices. 

In three dimensions the analogous problem is more difficult since the matrix 
no longer has independent random entries. It belongs to this new class of random 
matrices alluded above, for which large size asymptotics are not known. In this 
paper we initiate the study of such matrices by large deviation techniques which 
ultimately rely on Feynman graphs estimates. We prove a theorem which states 
that the probability that such random matrices develop a large eigenvalue is 
sufficiently small, so that a single slice analysis again can be controlled. 

This single slice analysis can be considered either as a step in the program to 
analyze by the RG method the situation where the energy lies in the middle of 
the conduction band, or as a full result in its own right, but for an energy out of 
the conduction band and close to its border. In this last case our probabilistic 
bounds imply via a cluster expansion that the averaged Green's function has 
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fast decay (for a Cfi° band cutoff, this decay is polynomial of arbitrarily high 
order), and that it is asymptotic to its perturbative expansion. 

3 Notations and Results 
3.1 Definitions 

We consider an electron in a weak random potential in three dimensions. The 
mass of the electron is set to m = 1/2 so that in these units the ordinary kinetic 
energy in Fourier space is simply p 2 . To model a conduction band, we consider a 
free Hamiltonian H which is simply a multiplication operator in Fourier space: 

Ho{p)=p 2 h(p) (1) 

where rj is some rotation invariant function which is 1 on a certain interval 
[a,b] and zero on an interval [a',b'}, with a' > 0. Therefore in this model Hq, 
the energy, is both infinite for large momenta \p\ > b' and for small momenta 
(\p\ < a'), so there is both an ultraviolet cutoff, as should be the case for 
any non-relativistic problem, and an infrared cutoff on small momenta. Lattice 
cutoffs could also be treated by our methods, replacing for instance the function 
P 2 by ~ 2 cos ^) (hence they are not rotation invariant). 

The interacting Hamiltonian is obtained by adding to H a random potential 

V: 

H = H + XV (2) 

We take V as a Gaussian ultra-local field in x-space, i.e. each V(x) is an 
independent Gaussian variable. More precisely, it means that V is distributed 
according to a normalized Gaussian measure d\i K of covariance k. We also 
include a volume cutoff: this means that we choose a large box A of volume |A| 
and side size L — IAI 1 / 3 , and restrict the operator V to that box. We can also 
use e.g. periodic boundary conditions for the operator H on the boundary dA 
(although this last condition is not essential) . In that case the Fourier-transform 
of a A-periodic function f(x) is 

f(p)= I e*-*f(x)dx (3) 

J A 

and verifies 

= Ta7 E e " ip - x /(p) with p e ( 2 -^f (4) 

We are interested in computing the averaged Green's function at energy E 
(with a Feynman prescription) as a kernel in x-spacc, or the density of states, 
which is the imaginary part of this kernel at coinciding points. Therefore we 
want to analyze 

< G± > = lim lim < G A . e > (5) 
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<G ^ > =j H Q + XV-(E + l s) d ^ ® 

In the following we will note < . > for integration with respect to V, and Prob 
for the probability of an event, measured with respect to d/i K . 

It is important to remark that we can put some kind of ultraviolet cutoff 
on the random potential V as well, because we are treating a non-relativistic 
problem. Usually this is done by considering the potential on a lattice, but we 
can also choose to stay in the continuum and to take the covariance k of V not 
to be a S distribution but a smooth positive type kernel with compact support 
or fast decay at the unit scale (the scale of our ultra-violet cutoff). 



3.2 Expected Results 

In the standard Anderson model at weak coupling in dimension 3 the first signif- 
icant open mathematical problem is to prove that for a conduction band of finite 
width and an energy E inside the band (that is E £ [a, b]) the thermodynamic 
limit of Green's functions defined above exist and has fast decay (exponential 
decay in the case of the lattice cutoff) for real small enough coupling A. One 
want also to prove that the decay rate behaves, as a function of the coupling, 
as: 

A = tt 2 A 2 +0(A 4 ) (7) 

All these conjectures stem from the well-known perturbative analysis of the 
model. It is based on expanding (6) as a power series in A (resolvent expansion), 
then integrating on V. This perturbative analysis is of course plagued by the 
usual divergence of perturbation theory. 

To obtain the leading behavior A ~ 7r 2 A 2 , it is sufficient to compute a self- 
consistent renormalization of the divergent tadpole graph, which leads to: 

A ~ A 2 9 f d d p . . . . (8) 

J V -(E + iA) w 

so that 

1 ~ a 2 f d d v ^ (9) 

J % 2 -£) 2 + A2 • W 

This equation gives the behavior (7) since E belongs to [a, b], where the cut- 
off r\ is 1. This self-consistent computation is the so-called "rainbow graphs" 
approximation. 

Thus the first main effect of V is to add a mass term in 0(A 2 ) which screens 

the singularity of the free propagator — — . Physically this means that the 

p z — E 

electron wave function loses memory of its phase by the incoherent reflections 
on V. 

To prove this fast decay of < G± > or cquivalently to analyze the behavior 
at small A of the density of states is our first natural goal in the small coupling 
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analysis of the Anderson model. Remark however that this result would still fail 
to distinguish between localization or the existence of extended states, which 
can be distinguished only by the decay properties of < G+G- >. Let us recall 
that in the localization regime (which in d = 3 means for A large) the function 
< G+(x, y)G- (x, y) > decays exponentially, whether in the diffusive or extended 
states regime (which in d = 3 should exist for A small), it is conjectured that 
the function < G + (x,y)G_(x,y) > decays only as \x — y]' 1 . To distinguish 
rigorously between these two regimes at d — 3, hence to establish the existence 
of the Mott- Anderson transition is the long term goal for our program. 

3.3 A Result 

Here we do not even prove yet this fast decay of < G± > when E is inside the 
band, but establish a weaker result. We consider an energy E not inside the 
band but near its edge. Also the coupling constant A is taken real and small 
enough so that both the width of the band b' — a' and the distance between 
the energy and the border of the band (say, a' — E for an energy close to the 
lower edge of the band) are large compared to A 2 . Under these conditions we 
shall prove that the averaged Green's function has fast decay on a spatial scale 
of order r _1 , where r = inf{6' — a', a' — E}. More precisely we can state the: 

Theorem I Let r = mf{b' — a', a' — E}. There exists some small constant e 
such that if X 2 < er, then < G± > exist (the convergence of its kernel as A — > oo 
and e — > being uniform on all compact sets) and for all integer m there exists 
a constant K m such that: 

I < G ± > (x, y)\ < K m —— — p — . (10) 

(1 + r\x — y\) m 

Furthermore this averaged Green's function is asymptotic to all orders its pertur- 
bative expansion in A. This expansion is obtained by expanding the resolvent (see 
(19) below) and performing the Gaussian integration onV in < G± > (x,y). 

This result in d = 2 is essentially the subject of [P]. The purpose of these 
notes is to extend it to d = 3, and explain the difference. The proof combines 
two aspects, a large/small field probabilistic analysis which is given here in 
detail, and a cluster expansion which is a rather standard method ([R][P]), so 
not repeated here. 

3.4 Next steps... 

The same result can be proved to hold also for a band width no longer small 
but finite, of order 1, the distance of the energy to the border of the band being 
still small (hence A 2 < er = a' — E); since the proof is technically slightly more 
difficult, we do not include it here. 
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Finally the much more difficult problem of an energy within the band, is 
somewhat analogous to spontaneous mass generation in field theory or in the 
BCS theory. Our result which concerns the energy outside the band is never- 
theless a relevant step towards this more difficult result, since one can expect to 
analyze the case of the energy inside the band by cutting slices closer and closer 
to it, as in a renormalization group analysis [P]-[FT]. We must however add 
that the situation is quite different in dimensions 2 and 3. In two dimensions 
we have a precise scenario in progress to solve this problem of an energy within 
the band, using the specific two dimensional rules of conservation of momentum 
[MPR]. In three dimensions we have not yet such a precise scenario on how to 
attack the problem of an energy within the band. 



3.5 Complex translation and limit e — > 

We can make an a priori partial renormalization and generate a non vanish- 
ing imaginary part in the denominator thanks to a complex translation of the 
integration contour of V, using the formulas 



J 'dfi c (X)T(X) = e 



{aCa)/2 



diic(X)T(X + iCa)e- iaJi . (11) 

These formulas are easily proved for T polynomial and then extended with a 
density argument to meromorphic functions provided we meet no pole. 

5 

Let us translate V by — i — 1, and get: 

n / „i51.V/\ \ 

^- W,/a L- {E + i s + ic) + >v ) < 12 > 

(The exact value of 5 will be fixed later on). Now we can take the limit e —* 
which yields 

/ c iSl.V/X \ 

°'-^< TOw )' (13) 

On the lattice this formula is exact but in the continuum it has to be slightly 
corrected to take into account the ultraviolet cutoff on V, but this is inessential 

n 

The cost of this operation is the bad factor e' A l 5 l 2X that we may have to 
compensate if some absolute value in the bounds cancels the oscillating term 
e iSi.v/\_ However let us explain the strategy: space will be divided into so-called 
large/small field regions. In the small field regions, perturbation theory works. 
The eigenvalues of XV are too small to cancel H — E. Therefore the imaginary 
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term iS is not really needed and the complex translation can be reversed, so that 
the bad factor has not to be compensated. In the large field region, this is no 
longer true, but the bad factor is more than compensated by the probabilistic 
factor associated to the measure dfi K which is very small. 



3.6 The resolvent expansion 

From now on the number r is chosen of the form r = M~i , where M is some 
fixed number and j a large positive integer. These are the usual notations which 
come from the renormalization group point of view [FT] , [FMRT] [P] . We also set 
E = 1 for simplicity. Let us also define rj(p) = [r]j(p)] 3 (the subscript j simply 
recalls that the width of the support of the function r\ is M - - 7 ). We introduce 
the free Green's function 

T={H Q -{l + i8))- 1 = Vo C Vj ; (14) 

where ^ 

C7= P 2 -^(l+i<5) • (15) 
Then we can factorize C and rewrite G as 



r _ |A|6 2 /2A 2 / m.v/x 1 c 



= A(e i61 - v ' x n{V)^ (16) 
where we used the operator product conventions 

(AB)(x,y) = j A(x,z)B(z,y)dz 
^ (AB)(p,q) = ±J2A(p,k)B(-k,q) (17) 

where 

A(p,q) = je vpx - iqy A(x,y)dxdy. (18) 

Because of its smooth covariance V is almost surely bounded. Thus for a 
given V, "R-{V) is analytic in A in a small domain around the real axis. This 
implies that 1Z is the unique analytic continuation of the operator series 

00 

S = Y,(-\) n m{ CHC...HC) Vj (19) 

™=° nH's 

with H = rjjVrij . (20) 
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4 Phase space Analysis 



We want to study G by means of a phase space analysis. In order to do so, 
we will cut the propagator C according to cells (called sectors) in the space 
of angular coordinates of the momentum p/\p\ (see subsection below). These 
sectors are also the main tool in the constructive renormalization group analysis 
of [FMRT1-4] . Space will also be cut according to a lattice of cubes of size 
M- 7 = R~ x , which is adapted to the decay of the propagator, or of the cutoff 
r]j. In momentum space V, seen as an operator, is V(p, q) = V(p — q). 



d=2 



d=3 




Figure 1 

The key difference between the constructive or phase space analysis in 2 and 
3 dimensions is that in two dimensions when p and q have the same fixed norm 
(here close to E = 1, since \p\ — 1 < 2r = 2M"- 7 ), the sump+q defines in a unique 
way the pair {p, q}, whether in three dimensions it defines it only up to a rotation 
by an angle (see Figure 1). The operator V can be considered as a matrix 
whose entries are labeled by this discrete set of sector indices. In this point of 
view, the geometrical behavior depicted in Figure 1 translates essentially in the 
following property: the matrix obtained in d = 2 is approximately an orthogonal 
ensemble with Wigner's semi-circle law as asymptotic behavior of the eigenvalues 
density [P], whether in three dimensions the matrix is of a constrained type, 
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and the asymptotic distribution for the eigenvalues can be probed only more 
indirectly, for instance, as will be shown below, through large order deviations 
which amount to the computation of particular Feynman graphs. We go on to 
make now this statement more precise. 

4.1 The discretized "single cube problem" as a Random 
Constrained Matrix 

The usual point of view in localization theory is to treat space as fundamental 
and to expand the off diagonal part of C, which couples different points. But 
for diffusion theory and extended states regimes, it is better to keep momentum 
space as fundamental. Since C is a diagonal operator in momentum space, the 
important part of the problem is to understand the behavior of the random 
operator H in momentum space. We show now how to discretize this operator 
as a matrix to understand better its qualitative behavior. 

From now on we restrict ourselves to dimension d = 3. Suppose that we 
cut the momentum shell support of the cutoff rjj into cells a (isotropic sectors) 
which have size M~ 3 in all directions: r]j = ^^Vj,^- There will be roughly 
N = M 2j such sectors (neglecting factors such as Air). Consider the covariance 
k of V. Since V is almost ultralocal, this covariance is almost 1 on a ball of 
large radius. But by momentum conservation, the operator H ajJ i = Tjj i(T Vrij i(T i 
is zero unless the momentum of V , which should be approximately a — a' , lies 
inside the ball of radius 2 + 4r, so approximately the ball of radius 2. In this 
ball k(j>) = 1 with a very good approximation. Cut the interior of this three 
dimensional ball of radius 2 into about M 3 ° ~ N 3 / 2 cells c of side size M~ J : 
k = k c . There will be an associated decomposition of V as a sum of random 
variables V c , each independent and distributed with covariance n c . 

When this decomposition of momentum space both for the electron cutoff 
rjj and for V has been made with smooth functions, there is fast dual 
spatial decay of the propagators rjj <a and k c on length scales M 3 . Therefore 
the heart of the theory is given by the problem restricted to a single spatial 
cube of side size M 3 . To understand this single cube theory we can neglect now 
the space dependence of rjj !tT and k c , keeping only the constraint of momentum 
conservation. 

In this approximation the different random variables V c (which are no longer 
fields but real random variables) become identically distributed. The matrix 
H has now N 2 entries, namely all pairs (a, a') where a and a' are isotropic 
sectors. We understand that since this N by N matrix depends on N 3 ^ 2 in- 
dependent variables only, it cannot be of the free type usually considered in 
theoretical physics (orthogonal, hcrmitian or symplectic ensemble). It is of a 
new constrained type, and the constraints stem from the geometry of momen- 
tum addition on the singular manifold (here the sphere). These constraints are 
related to the absence, in dimension 3 or more, of conservation rules for diffusion 
and of completely integrable models. Because of these constraints, the matrices 
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obtained are not invariant under a group such as the orthogonal, unitary and 
symplectic group; the statistics do no longer factorize and reduce to eigenvalue 
statistics (namely there are "preferred eigenvectors" as well), orthogonal poly- 
nomial methods do not apply, etc., and nothing seems to be known up to now 
about such matrices. 

These matrices can nevertheless be described and studied numerically on 
a computer very easily with the following program (here we simplified again 
slightly the initial problem, which involved complex matrices, replacing it by an 
analogous real problem with sums of momenta instead of differences): 

Take a small number r = M~i . Divide the unit sphere into TV ~ M 2j 
sectors of approximately similar area 47rr 2 . Divide the ball of radius two into 
TV 3 / 2 ~ M 3j / 2 cubic cells of volume approximately (4/3)87rr 3 . Generate with a 
random number generator TV 3 / 2 independent random variables V c , one for each 
cell of the ball. Associate to any pair of sectors (a, a') the unique cell c = c(er, a') 
such that the sum of the two momenta at the center of a and a' fall into the 
cell c, and plug the corresponding variable V c (a,a') at the corresponding entry 
of the matrix, that is consider the (real symmetric) N by N matrix H(a, a') = 
V c (a, a'). Diagonalize it, plot the eigenvalues, and repeat the operation a lot of 
times to get some statistics. 

We have implemented this program and run it on a computer for N up 
to about 600 [G. Poirot, unpublished]. We cannot yet predict clearly scaling 
laws or tail estimates from this numerical study, but the results indicate that 
the eigenvalue distribution is differently peaked and extends more towards large 
eigenvalues that in the free case. This is not surprising, since constraints typi- 
cally increase both the degeneracy hence the small eigenvalues of a matrix, and 
allow for larger exceptional eigenvalues. For instance the completely constrained 
matrix with the same random value at every entry has N — 1 eigenvalues 0, and 
a single large eigenvalue of size N. 

We close this parenthesis on this simple discretization of the random matrix 
model, and return now to the true problem which is complicated by spatial de- 
pendence. For this true problem one gets better bounds by considering angular 
sectors, called anisotropic, which are wider in the tangential than in the radial 
directions [FMRT1]. We introduce now these sectors and the precise tools to 
prove Theorem I below. 

4.2 Anisotropic Sector Decomposition 

We noticed that the cutoff r]j implies the condition: 

||p|-l|<2M- j . (21) 

We construct a partition of the unit sphere into anisotropic sectors [FMRT2] 
S a , by projecting e.g. a standard division of the faces of the cube [— 1, l] 3 into 
6 x 4M( 2j ')/ 2 = 24AP square plaquettes of side M~ j / 2 , hence area M~ j , onto 
the sphere. This gives a sharp partition of the sphere into 24AP "spherical 
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plaquettes" or sectors a of variable area (between M i and , since the 
maximal distance of the cube to the origin is y/3). We call k a the center of 
mass of a. From this sharp partition one can build an associated smooth C°° 
scaled partition of unity \a on the sphere with Urisohn's lemma. Finally we 
transport this partition into a partition of the support of the cutoff r]j by writing: 

»fe=I>a (22) 

a 

where r] a (p) = r]j(p) ■ Xa(p/\p\)- 

The support S a of a sector a, or of an rj a function, is therefore included 
in a parallelepipedic slab, of center k a , which has small thickness I = M~ J in 
the direction parallel to k a and which is longer in the two other perpendicular 
directions where it has length L = 4M^'/ 2 (1 + jj) = 0(1)M~^ 2 . 

The cube and our sector partition is invariant under parity: p — » —p. There- 
fore for any sector a with center of mass k a we can define the antipodal sector 
a = —a. All other sectors have support at distance at least say (l/4)M~-7/ 2 of 
the antipod center k & = —k a , if the partition of unity \a does not spread too 
much onto the nearest neighbors, a fact we can assume from now on. 

Then if (3 ^ a, we define 

k a -(3 = k a -k p = 27 Q _ /3 cosx Qj/3 (23) 

where x a .p — 9 a ,/3/2 is half the angle 6 a j) (modulo 7r) between k a and — hp, so 
x a ,/3 € [0, 7r/2[, and 7 a -/3 is a well defined unit vector. 



4.3 Triplet of cubes operator 

It is natural to consider a smooth partition of the space by a lattice V of cubes 
A of side M- 7 , with "characteristic functions" \a ■ There is a corresponding 
orthogonal decomposition of the random variable V = V A , by writing its 
covariance n(x, y) — J2 A k 1 / 2 (x, z)xa(z)k 1 ^ 2 (z, y). 

We define H A * (x, y) = T]jV A *rij and for each triplet A = (A, A', A") e D 3 
or each pair (A, A') e V 3 

H a( x 'V) = H A " = Xa(x)H a " (x,y)xA'(y) 

A A' 

+ XA'(x)H A "(x,y)xA(y) (24) 

Similarly we define 

C(A,A')(x,y) = XA(x)C(x,y) X A'(y) 

+ XA>(x)C(x,y)xA(y) (25) 
to keep the operators symmetric. 
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5 Size of the operator H 



Let us repeat that because we are in dimension d — 3, the operator H localized 
in a box, considered as a matrix whose entries are the sectors of ingoing and 
outgoing momenta, does not, even approximately, look like a classical ensem- 
ble such as the orthogonal or Hcrmitian ensembles of random matrices which 
appeared in dimension 2 [P] . So Wigner's semi-circle law and the usual asymp- 
totics do not apply to this problem. We must use a different, less accurate, 
method to bound the tail of the distribution, namely the probability that some 
eigenvalue becomes large. 

Let Nj = M 2 i, and 

D(A) = (1 + M~ j d(A, A') + M~ j d(A', A")) (26) 
Our main probabilistic bound is 
Theorem II 

Let p > 1 be some integer. There exists a constant C p > 1 such that for 
A>C P 

Prob{||ff A || > A.DiAyPN- 1 /*} 
< (AD(A) 2p )~ N3/16 . (27) 



The proof is quite lengthy, and is the subject of the rest of the paper. 

Recall that for a compact symmetric trace class operator such as H (it is 
easily proved trace class because of the cutoffs), the ordinary norm is bounded 
by {trH 2n ) 1 / 2n for any positive integer n. Therefore 

Prob{||ff A || > A.DiAyPN- 1 ^} 

< (AD(A)-JW- 1 / 4 )- 2 ™ J d^(V A ")trH 2n . (28) 

Performing the Gaussian integration we obtain a sum over all Feynman 
graphs obtained by contracting 2n insertions of V A on a loop of 2n propagators 
which are r]j(XA + XA') 2 1j- n 1S some integer that will be adjusted later. 
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Figure 2: A graph G 
and its associated cf) 4 graph 
obtained by contracting the dotted lines 



More precisely 

f d^(V A ")trH 2n = £ I G , (29) 
Geg m 

where the amplitude for the graph G is noted Iq. 

To evaluate the amplitude of a given graph Ia,we split the spatial integrals 
into cubes of side size M^ 2 , hence volume iV 3 / 4 . In the language of the </> 4 field 
theory, the "vertices" correspond to the contractions of the initial V random 
field. They are propagators W(x,y) = k 1 / 2 (x, z)xa(z)k 1 ^ 2 (z, y). They are 
represented by the dotted lines in Figure 2. There are also the "ordinary lines" 
which correspond to the electron propagators rjjXAvA'Vj- They are represented 
by solid lines in Figure 2. 
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We split these electron propagators as Ta,A' = J2 a a' ^AA'VaXAuA'Va- 
Since the functions xa are scaled and smooth, we have momentum conser- 
vation, which means that a must be close to a', and then we have spatial decay 
of r dual to the length of the cell. Therefore: 

Lemma 1 

For any integer q, there exists a constant C q such that 

C q M- 2 3 



< 



[1 + M-^d 2 (x, A U A')] [1 + Mid 2 (a, a')] 
1 



1 + d 2 (x, y) [1 + M-i d 2 a± (x, y)] Q 



(30) 



where d 2 (a,a') is the distance on the sphere between sectors a and a', and d/i 
and d± are the distances in the axis parallel to the center of a and in the plane 
perpendicular to it. 

Proof This is an easy exercise in integration by parts, similar to Lemma 4 in 
[FMRT2] (pg 689). 

Applying also standard rules of integration by parts (see e.g. [FMRT2]), we 
can extract a very small factor if sectors a\, «2, «3 and a 4 meet at a given 
vertex, and do not contain points which add up to 0. This mean that we can 
bound every vertex, extracting such a momentum conservation factor, plus some 
distance factor, by 



[l + A'P|di+a 2 +a3 + a4|] 9 
1 



[l + M-^d 2 (z„,A")]' 



(31) 



where on is the center of the sector on . 

Substituting these bounds we obtain our general bound for a graph of order 
n = 2m, taking into account the fact that each half vertex joins successively A, 
A" and A' or A', A" and A: 

Iq < cl n D{K)- qn 



2_] / dxi / dx 2 ---dx n 

Tfr,-, JAVA' Jr 3 



{ai.leG} 

n 



[1 + M3\a 1 +a 2 + a 3 + a 4 \ 
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n 



N~ 



\ [l + M^ 2 {a u a' l )\ 



l + M-jd 2 ai± ( Xl ,yi) 



(32) 



where the sum runs over attributions of sectors to the lines, and spatial inte- 
grations over the vertices have been performed, except at one point per vertex. 

Now by a combinatoric factor 2 per vertex we can decide if there exists at 
each vertex a sum or a difference between two incident sectors which is small in 
the sense that \a v ± a' v \ < M - - 7 ^ 1 / 2 ) -0 ). a is some small number that we will 
fix later to a = 1/32. In this way we cut the amplitude Iq in 2" elements as 
Iq = ^2t^G,t, where the sum over t runs over the choices at each vertex whether 
there is or not such a small sum or difference. 

Vertices such that there is no such small sum or difference are called "twisted" 
or non planar. The vertices with such a small difference are called "untwisted" 
since their four momenta typically will lie approximately in a common plane, or 
otherwise the vertex factor in (31) will be very small. The number of twisted 
vertices is called n u and the number of untwisted vertices is n u , with n = n t +n u . 

Now for each "untwisted vertex" , we choose the particular pair of momenta 
realizing the smallest sum or difference (or an arbitrary one if several pair realize 
this minimum) and we split the vertex into two half vertices and an intermediate 
wavy line according to this splitting. The graph has then n t quartic vertices 
and n u pairs of trilinear vertices. Then we cut all the corresponding wavy lines. 
The graph splits then into c "ordinary" connected components which contain 
at least one twisted vertex, plus possibly other twisted vertices and chains of 
untwisted half vertices, and into d pure "cycles" of untwisted half vertices, with 
c + c' > 1. The result of this splitting process is pictured in Figure 3, for a graph 
with 15 vertices, 8 of which are untwisted and seven twisted. It splits into c = 5 
connected components if the wavy lines (untwisted vertices) are removed. We 
have c' = 3 of these 5 components which are cycles which do not contain twisted 
vertices. 

We select a tree of wavy lines relating together all these components, and will 
perform the spatial integration and attribution of momentum sectors according 
to what we call "painting rules" for each connected components. The tree of 
wavy lines is ordered by chosing a particular root, which correspond to a special 
connected component G (either ordinary or cycle). 
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<\A/W Untwisted vertices 
• Ordinary (twisted) vertices 



Figure 3: Splitting a graph 
into twisted and untwisted vertices 

The standard weight for integration of a vertex with the help of a single 
propagator hooked to it is the volume of a tube of size M^ 2 in the two per- 
pendicular directions to a sector, and M- 7 in the parallel direction, hence it is 
N. However if the vertex is twisted and integrated by taking into account the 
decay of two propagators hooked to it, remark that the corresponding weight 
is at most iV 1-a / 2 , since we gain the fact that the integration is located at the 
intersection of two tubes who make an ang le at least M"^ 1 ^)-"). Finally 
there is one single vertex in the root component whose integration costs the full 
volume of A U A U A", hence M 3j = TV 3 / 2 . 

Now let us perform momentum attributions and vertex integrations. We 
perform these attributions in an order induced from the tree chosen and from 
fixing rules associated to each connected components. Let Gk be pure cycles 
of pk untwisted half vertices. Let r*, be the number of such half-vertices whose 
other half vertex is not yet fixed in the process (i.e. is lower in the tree). We 
have rfe > 1, except possibly for the last cycle G if it is the root of the tree. 
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Then the first sector attribution to a line of Gk is made on a half vertex of the 
type rk (except possibly if k = 0). It costs TV 1 / 2 (total number of sectors). Each 
next attribution costs cither N a (number of sectors in the cone of opening angle 
_M-j(( 1 / 2 )-a)) ; or 0(1), depending whether the sectors of the other half vertex 
to which this half vertex is associated has been already attributed or not. So the 
total factor is jV^+afo-i) = N^/2)-a+ar k ^ except possibly for k = 0, where 
it can be N 1 / 2 if r = 0, so there is an extra N a in that case, with respect to 
the formula jV( 1 / 2 )-° +or *. 

Finally it remains to attribute the momenta in the ordinary components. 
Let us contract the chains of untwisted vertices to single lines (but with an 
index p for the number of insertions of half untwisted vertices, and an index 
r for the number of such insertions whose other half vertex is not yet fixed). 
Fixing all the sectors of such a line is like fixing one of them at the end, plus a 
factor N a - r . 

To complete the bound, it remains therefore to compute the weight of fixing 
sectors in a vacuum connected <fi 4 graph solely made of twisted vertices. Let G 
be such a graph, and s its number of vertices. 

A fixing rule P on G is an ordering of the vertices of G as v\, v s such that 
for every k > 1 there exists a line hooked at one end to Vk and at the other 
end to one vertex among {v\, ...,Vk-i}- There always obviously exist such fixing 
rules. For instance to any tree T of G and any way of "turning around the tree" 
(i.e. choice of a root and relative ordering of the non trivial forks going up in the 
tree) is associated a fixing rule corresponding to the ordering in which vertices 
of G are met by "climbing around the tree" in that way. 

The fc-th vertex Vk for the fixing rule P is called of type 4, 3, 2, 1 or with 
respect to P if the number of half lines hooked at one end to Vk and at the other 
end to a vertex not among {v\, Vk-i} is 4,3,2,1, or 0. The number of vertices 
of type r is called n r . Obviously + + 712 + n\ + no = s, and = 1 since 
the first vertex of the fixing rule is the only one of type 4. 

To perform the sum over the values of sectors in the graph we fix first the 
indices of vertex V\, which is the only one of type 4. This costs TV 1 / 2 ^ 1 / 2 ^ 1 / 4 = 
iV 5 / 4 . Then vertices of type 3 cost V 1 / 2 ^ 1 / 4 = N 3 ^. Vertices of type 2 cost 
N 1 ' 4 , because they are twisted. Vertices of type 1 or zero cost nothing. Finally 
we should take into account the fact that because of the twisting, vertices of 
type 2,1 or 0, when integrated, are integrated with the decay of two propagators 
whose sectors are known and have a minimum angle. This allows a gain of N~ a , 
as explained above, for the spatial integration of these vertices. 

Collecting all factors, the key point is to remark that the factors N arh and 
N ar combine into N ariu , since one of the two half vertices of any untwisted 
vertex is always fixed before the other. Adding the extra N a eventual factor 
for Go, we arrive at the following estimate, where we recall that n u is the total 
number of untwisted vertices: 

let < C 3n D(A)~ qn N~ 2n N n+1 / 2 N a< - n < 1+1 ' 1 
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jy(l/2-o)c' + (5/4)c+(3/4)n 3 +((l/4)-o)n 2 -a(ni+n ) ^33) 

where we have n = n u + c+ns + n2 + ni+no, and 2n — 2n u + 4c+3n3 + 2n2 + ni, 
counting respectively the vertices and the lines. Taking into account the second 
relation we find 

let < Cg n D(A)~ qn N 1 < /2 ~ n,/2 N c> V 2 -™../ 2 

jyc/4-ra 2 /4-ni/ 4 jya(n u +l-(n +ni+n 2 +c')) (34) 

Taking q = 4p, the factor 7V-™/ 2 £>(A)- W 2 kills the N n ' 2 D{A) 2 P n in for- 
mula (27). It remains: 

N n / 2 D(AfP n \I G , t \=J git 
< D(A)~ 2pn N a+1 ^ 2 N c 'I 2 

^y— ra„/2+c/4— n 2 /4— ni/4jy<m„ — a(na+ni+n2 + c') (35) 

From the fact that there is a tree of untwisted "wavy lines" between the con- 
nected components of type c and c', we know that c + c' < 1 + n u . Therefore 

J g ,t < K^D{K)- 2pn N a+:i ' A 

jY-(l/4-o)(n„-c')-n 2 /4-ni/4jy-a(rto+ni+n2) (36) 

Choosing a=l/16, we have: 

J g ,t < K%D(A)- 2pn N a + 3 / 4 

jy- (3/16)(n„-c')-(l/16)(n +ni+n 2 ) (37) 

We have n t — c + n 3 + n 2 +«i +n , and 2n t = 4c + 3n 3 + 2n 2 + ni, therefore 
2c + n 3 = m + 2n . Therefore c + n 3 < m + 2n , so n t < n 2 + 2m + 3n < 
3(n 2 + ni + no). Hence 

J g ,t < K™D{A)~ 2pn N a+i / 4 

jy-(3/16)[(n u -c')-(3/16)n t ] 

< i^ l Z)(A)~ 2p "N a+3 / 4 iV _(3 / 16)( ™ _c ') (38) 

Now we consider the following reduction process. We contract a certain 
number x < c! of pure cycles in the following inductive way: first we contract 
every tadpole (pure cycle with p = 1), i.e. suppress it together with one of 
the two external lines hooked to the tadpole. Then we repeat this operation 
inductively (since the contraction at first stage can create new tadpoles), until 
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complete exhaustion. In the case of the graph of Figure 3, there is a single 
tadpole so the operation stops after one step, but in general it can be more 
complicated (see Figure 4). 




Figure 4: The contraction of tadpoles 




Figure 5: a pure loop of 8 bubbles 



Let x\ be the number of cycles contracted in this way. This contraction 
process continues until either we arrive at the empty graph because at some 
stage we obtained a graph with c' — 2 made of two tadpoles, which we then 
delete or until we arrive at a graph with d — x\ cycles without any tadpole. 
In this last case either the graph is a pure loop of x 2 bubbles that we destroy 
(see Figure 5), or it is not. In this second case, the graph may still contain 
pure bubbles, i.e. loops without twisted vertices of length 2. For instance in the 
case of the graph of Figure 3, after contraction of the tadpole, there is one such 
bubble which appears. 
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We perform then a second operation of contraction, which contracts every 
maximal chain of such bubbles (pure cycles of length 2 without twisted vertices) 
to a single new type of wavy line. This second operation has not to be repeated 
inductively in contrast with the first, since it cannot regenerate any bubble (see 
Figure 6). (Wavy lines have been omitted in Figure 4-5-6 for simpler pictures). 




Figure 6: the reduction of bubbles 



Let x 2 be the total number of bubbles contracted in this way. We arrive at 
a graph with d — x\ — x 2 cycles, each of them having at least three lines. 
Let x = x\ + x 2 . 

Lemma 2 We have x = x\ + x 2 < /2(G) + 2, where /2(G) is the cardinal of 
the largest forests of connected 4-point subgraphs of G, defined in [CRJ. 

Proof: Obviously each contracted cycle can be associated to a new 4 point 
subgraph in a forest, except if a; = c', in which case when one destroys the 
two tadpoles there is no new 4-point subgraph associated, or when one destroys 
the loop of k bubbles, there is only k — 1 corresponding 4-point subgraphs in a 
forest. This two cases are exclusive, which explains the 2 in the lemma. 

We have obviously x < c' hence n — c' > n — x > n — f 2 (G) — 2. We use now: 

Lemma 3 (particular case of [CR], theorem II) The number of Wick 
contractions of half-vertices giving rise to a vacuum connected graph with n 
vertices and f 2 {G)+2 = f, is bounded by K n n\/ /! < K n n n ~f , for some constant 
K. 

Proof See ([CR], Theorem II). 

Combining this, with the previous bound we obtain (since / < n + 2 ([CR])): 

G 

J2 {KK q ) n n n - f D{A)- 2pn N a+3 ^N~ < - 3 ^ 16 ^^ 

f<n+2 
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< C n D(A)- 2pn N a+3/i (n/N 3/w ) n - f (39) 
where C = 2K.K q is a constant. We choose now n — TV 3 / 16 and obtain: 

Prob{||# A || > A.D(A)- p N- 1/4 } 

< A- 2nVW D(A)- 2 p n3/W N -+3/*c n3/16 (40) 

If A 1 / 2 is larger than C, we obtain easily the bound (32) of the theorem 
(since A~ n3/W / 2 7V a+3 / 4 < 1, A being large and N > 1, and we can also choose 
A such that CA^ 1 / 2 < 1). We have therefore in that case: 

Prob{||ff A || > A.D(A)~ P N~ 1 / 4 } 

< A- N3/16 D(A)- 2pN3/16 (41) 
which completes the proof of Theorem II. ■ 



In order to obtain a theorem no longer on H but on the full operator XCH 
that occurs in our resolvent expansions, we need to add to the estimate of 
Theorem I the power counting coming from the coupling constant and the norm 
of the propagator. This adds a factor eN 1 / 4 to the norm of H. The factor e 
comes from the distance to the last slice, since recall that in section 3 we assumed 
that A 2 < eM -1 . Finally we obtain, taking e = A~ 2 small, the following 
estimate: 

Theorem III 

// the coupling is kept small enough (i.e. e < A~ 2 , where A is defined by 
Theorem II) we have: 

Prob{||ACH A || > A^DiA)-?} 
< A- N3/16 D(A)- 2pN3/ie (42) 

To complete the proof of Theorem I, one combines the estimate of Theorem 
III with a large versus small field cluster expansion. This expansion is exactly 
similar to the one of [P], to which we refer the reader. The decay factor D(A)~ P 
is roughly speaking needed to sum the various triplets A with respect to one of 
their elements, and the factor A~ N ^' is sufficient to compensate the bad factors 
due to the imaginary translation (see (13)) which come from the estimation of 
the resolvent in the large field regions. 
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